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Abstract: We investigate in detail a model that contains an additional SU{2) singlet and 
triplet scalar fields than the Standard Model (SM). This allows the radiative generation 
of Majorana neutrino masses at two-loop order with the help of doubly charged Higgs 
bosons that arise from the extended Higgs sector. The phenomenology of the Higgs and 
neutrino sectors of the model is studied. We give the analytical form of the masses of scalar 
and pseudoscalar bosons and their mixings, and the structure of the active neutrino mass 
matrix. It is found that the model accommodates only normal neutrino mass hierarchy, 
and that there is a large parameter space where the doubly charged Higgs can be observed 
at the Large Hadron Collider (LHC), thereby making it testable at the LHC. Furthermore, 
the neutrino-less double beta decays arise predominantly from exchange processes 

involving the doubly charged Higgs, whose existence is thus unmistakable if ^v(3(3 decays 
are observed. The production and decays of the doubly charged Higgs are analyzed, and 
distinct and distinguishing signals are discussed. 
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1. Introduction 

The origin of small active neutrino masses remains one of the most challenging problem in 
physics. The small neutrino masses generated through the seesaw mechanism is popularly 
viewed as heralding new physics at scales larger than 10^^ GeV, and thus provide a win- 
dow to Grand Unified Theories (GUTs) with or without supersymmetry. Crucial to the 
construction is the introduction of heavy Standard Model (SM) singlet fermions commonly 
known as sterile neutrinos. 

Recently, the idea of extra spatial dimensions together with brane world scenarios offers 
a very different perspective to the question of neutrino masses. Here their smallness results 
from either the suppression factors associated with the relatively large extra dimensions, 
or from the small overlap between the wave functions of the sterile neutrinos in the extra 
dimensions. 

It is interesting to note that these different perspectives can be incorporated into a 
single framework in brane world scenarios; a recent discussion can be found in Q. However 
the existence of sterile neutrinos is required in both constructions which, along with the 
value of their masses, are all important questions by themselves. To date the best informa- 
tion on light sterile neutrinos comes from cosmological considerations; direct experimental 
tests are very challenging due to the fact that they have no SM interactions. 

It is well known that the masses of active neutrinos can be generated without sterile 
right-handed (RH) neutrinos via quantum loop effects. Without the RH states there are 
no Dirac couplings of the SM lepton doublet to the Higgs fields, and consequently the 
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active neutrinos can only have Majorana masses. The prototype model was constructed 
in 1^] where there is an extended Higgs sector, and the gauge symmetry is that of the SM. 
Crucial to the construction was the use of an SU{2) singlet Higgs field with a nontrivial 
hypercharge. Unfortunately the model gives rise to bimaximal neutrino mixings which is 
disfavored by the most recent neutrino data (for a recent review see |^]). More realistic 
neutrino masses can be obtained using doubly charged Higgs fields 

In our construction, we keep the SM gauge group and extend the Higgs sector by 
adding both an SU (2) triplet and a doubly charged singlet field. We also postulate that 
lepton number violating effects take place only in the scalar potential, while the rest of the 
Lagrangian respect lepton number. A brief discussion of our model has already appeared 
in where we showed how naturally small neutrino masses can arise from just two- 
loop radiative corrections. In this paper we will give a detail discussion of rich scalar 
phenomenology of the model. In particular the signals at the Large Hadron Collider (LHC) 
are investigated. 

The rest of the paper is organized as follow. In Sec. 2 we describe in detail our model. 
We work out the constraints on the vacuum expectation value and other parameters that 
control the mass of the two physical doubly charged Higgs bosons Pi^2^- We show that 
at least one of the doubly charged Higgs can have mass at the electroweak scale if we 
demand the theory be perturbative up to the TeV scale. In Sec. 3 we discuss in detail the 
phenomenology of the neutrino sector in our model. We examine closely the neutrino mass 
matrix and the constraints from the oscillation data. We show that normal hierarchy arises 
naturally in our model, and we place constraints on the neutrino-lepton Yukawa couplings. 
Lastly, we discuss the implications these constraints have on the Qv(3(3 decays of nuclei in 
our model. In Sec. 4 we discuss the phenomenology of the doubly charged Higgs production 
at the LHC, and their decays. We show that the decay pattern of the doubly charged Higgs 
in our model can be very different, and can therefore be used to distinguish our model from 
others that also contain doubly charged Higgs. Sec. 5 contains our conclusions. 



2. A minimal model with radiative neutrino mass generation. 

The model is based on the SM gauge symmetry with an extended Higgs sector and minimal 
matter content. Group theory dictates that only SU{2) singlets and triplets are allowed 
for the generation of Majorana masses for the neutrinos. Besides the SM Higgs doublet 
given by 




:ir20, (2.1) 
we introduce a complex triplet Higgs T represented by a 2 x 2 matrix 

(2.2) 



^ V2 



as well as a complex singlet scalar The subscripts denote the weak hypercharges of 

the fields as given by the relation Q = + ^Y. The most general potential for the scalar 
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fields is given by 

y((/.,r,^) = + A^(</.V)^ - f4Tr{T^T) + \T[Tr{T^T)f + X'TTr{T^TT^T) 

+ prr(ri'r^t^) + ^(^^r^i') - M{(j)'^T^(j)) + /i.c.j . (2.3) 

One can assign ^ a lepton number 2, and T a lepton number 0. Then terms in the square 
brackets at the end of Eq. ( |2.3D contain lepton number violating interactions. We will 
take both /xj. to be positive so that spontaneous symmetry breaking (SSB) takes place. 
Minimizing the potential gives us the VEVs: {(j)^) = and {T^) = 

With the additional fields T and ^, two new Yukawa terms can be constructed that 
are allowed by the gauge symmetry. The first is Yabl'^^j^lbR'^ i which is lepton number 
conserving. Here a, b are family indices and la is a lepton singlet. The second is LLT that 
violates lepton number, which we assume not to occur at the tree level *. The absence 
of this term frees one from having to put in by hand a very small value of vt in the eV 
range that plagues other Higgs models of neutrino masses. Adding in the SM terms and 
the covariant derivatives of T and ^ we have a complete renormalizable model. 



From Eq. (O) it can be seen that the various Higgs fields will mix among themselves. 
In particular, the pair DKecp^ and ^itT^ will mix that give rise to two physical neutral 
scalars, and P", and the pair 3m cj/' and JmT° will mix with one combination that is 
eaten by the Z boson to leave a physical pseudoscalar T^. Similarly for the charged states 
(j)^ and r^, one combination will be eaten by the W bosons leaving only a pair of singly 
charged scalars. Finally the weak eigenstates r=^=^ and will also mix to form 

physical states and P^^, with the mixing angle denoted by 6, and masses Mi and 
M2 respectively. All the masses and mixing angles are free parameters in our model, which 
we can use to replace some of the parameters in V{<j),T,'^). They are to be determined 
experimentally. In summary, the physical spin particles consist of a pair of singly charged 
bosons, P^, two pairs of doubly charged bosons P^^ and ^2"^=^, a pair of Higgs scalars 
and P^, and a pseudoscalar A^. 

The value of vt is constrained by the electroweak phenomenology. After the elec- 
troweak symmetry breaking, the W and Z bosons pick up masses at the tree level given 

by 

M^ = ^^{v'+2v't), M|= (v' + ^v't), (2.4) 

4 4 cos"' Uw 

where we have used standard notations, and the tree level relation e = gsm6\y holds. From 

the Particle Data Group (PDG) we have p = 1.002+;°°°^ |[ and Mw = 80.41 GeV §. 

This implies that vt < 4.41 GeV. We will see below that this is a controlling scale for the 

neutrino masses. 



'There are several ways to naturally suppress the Yukawa couplings of LaLtT. One way is to embed the 
model in a 5-dimensional set-up and compactify the extra dimensions on an orbifold Si /Z2 ■ The lepton and 
Higgs fields are then assigned with different orbifold parities to forbid the LLT term while still allows the 
lltp term. Another way is to further extend the Higgs sector by including a second Higgs doublet and then 
employ an appropriate discrete symmetry. The LLT term will be generated radiatively after symmetry 
breaking, but is small. 
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In the Umit where A, M — > the model conserves lepton number, and is thus technically 
natural. Since A is a dimensionless coupling, it is expected to be of order unity (< 4tt) 
so that perturbation is valid, which is assumed throughout this paper. The value of M is 
important in setting the scale of lepton number violation. It enters in the conditions for 
minimizing the scalar potential V{(j),T, 

-/i^ + X^v^ + ^K+i^T - V2MvT = , (2.5) 

-//| + A+4 + -/€+?;2_ _\ 0, (2.6) 



1 V2\VT 

where k+ = ki + K2 and A+ = Ay + A^^. Taking v to be of the electroweak scale and 
vt/v ~ 0.02, the interesting limiting cases are: 



A. M ~ vt- The minimum conditions, Eq. ( |2.5| ), can be naturally satisfied without fine 
tuning between the parameters if /Ut ~ ^t- 

B. M ~ f ^> Vt- To satisfy the minimum conditions, fjj^ ~ v'^ /vt is required. This 
appears to be unnatural although not forbidden. 

C. M > v: The minimum conditions can only be satisfied by tuning the dimensionful 
and/or the dimensionless couplings. We will not consider this case. 

For convenience define oj = . Then Case A and B correspond to w ~ 1 and 1 
respectively. Qualitatively, we see that a; < 1 is more natural in our model. We will 
therefore concentrate mostly in this region of the parameter space below. 

We now turn to the masses of the physical scalar and pseudoscalar particles in our 
model. The mass of the singly charged Higgs boson is given by 

M|,± = (a;-^)(t;2 + 24). (2.7) 

For cj ~ 1, we expect the charged Higgs to have mass in the 100 to 1000 GeV range. 
For the two doubly charged scalar s, their masses are given by 

1 



+ c T \/462 + (c - a)2 , (2.8) 



where the Pi state takes the upper sign, and 

1 11 

Note that m {m? > 0) is a mass parameter for the singlet and should not be confused with 
the physical mass. As such it is in general not constrained. 

Consider now Case A. In the limit where m is large (m ^ u), we have from Eqs. 
and (PD: 



^ AX + o(^), (2.10) 

4 



Ml^ = -{2UJ-K2)V 



Ml =m^ + + PVt) + O (%] . (2.11 
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We see in this limit, Mp^ saturates to an m-independent value, y^{2uj — K2)v'^ — X'rpV^, 
which is also its maximal value for a given set of model parameters. On the other hand, 
increases as m, which means the P2 state will be too heavy to be of interest to the 
LHC in the large m limit. 

Mp, (GeV) 




2 4 6 8 10 12 ' 

Figure 1: Maximum value of Mp^ for vt ^ M = 4 GeV, and |A^p| set to 47r. 




Figure 2: A/p^ as a function of m for |k2| ~ 0.5, 0.25, 0.125 in units of 47r, with vt = M = 4 GeV 
and A = — = 1. 



In Fig. H we plot the maximal value of Mp^ as a function of \k2\ = —K2, with vt and 
M set to 4 GeV. The coupling is set to |A^| = — A^ = 47r, the upper limit under which 
perturbation is expected to be valid. In Fig. |2| we plot Mp^ as a function of m for three 
different values of \n2\, with all dimensionless couplings kept perturbative. We see clearly 
here the saturation of Mp^ at large values of m. Figs. || and ^ show that the mass range 
of the Pi state is well within the reach of the LHC, its existence is thus a testable feature 
of our model at the LHC. 

Note that if we take the opposite limit where m ^ 0, two weak scale doubly charged 
scalars are possible. For example, if we take w = 1, A = 1, Kvjr = —K2 = 2, and p = 2A^ = 2, 
we get = 219.3 GeV and Mp^ = 365.8 GeV. However, this will not hold in Case B. In 
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this case the P2 state will be heavy with mass above a TeV, while the Pi state remains at 
the weak scale. 

The doubly charged scalars form a two-level system in which the mass and weak 
eigenstates are related by 



P. 



cos 5 sin 5 
— sin 5 cos 5 




(2.12) 



The mixing angle 5 is a measurable physical parameter. It is given by 

- {c-ar 



sin 25 



1 + 



462 



2Av2 



2A 



(2.13) 



where a, b and c are those given in Eq. (III). For Case A, if m2 < t;2, the mixing can be 
large and close to maximal. But if m?' ^ t;2, the mixing will be small, which is expected 
since the two states are widely split. For case B, large mixing can be achieved only if a 
cancellation occurs between the various parameter in Eq. (2.13). 

We now turn to the three physical neutral scalar bosons in our model. The pseudoscalar 
has mass given by 

M|o = ]:uj{v'^ + 44) . (2.14) 

Note that if w ^ 0, T° becomes a Major on. 

The masses of the neutral scalars hP and again have the general form 

1 r / 



hO,po 



+ c't \/46'2 + (c' - a')2 



where the state takes the upper sign, and 



b' 



-K+ -UJ ]VVT, 



c' = A+4 + 2*^^^ 



(2.15) 



(2.16) 



The physical neutral scalars also form a two-level system in which the mass and weak 



sin2i? 



1 + 



1 + 



angle. 


1? 


given by 


(c' - a' 




1 

2 


46'2 










2A+ 4 


164 v 




. - 2a; ?;2 



+ 



uj — 2Xa 



2uj 



(2.17) 



It can be seen that "d is of order vt/v for both Case A and B. 

For Case A, and P^ can both be light and almost degenerate. If they are lighter 
than half the Z boson mass they will contribute to its invisible width P] : 



r(z' 



6V27r 



. M2o+M2 y 
' Ml j 



(2.18) 
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where the notation is standard. Demanding that this contributes less than 150 MeV to the 
invisible width we obtain > 0.016. For Case B, all the neutral bosons have weak scale 
masses and the above limit does not apply. However, we still expect and to be close 
in mass. 

We summarize our findings on the masses of the scalar and pseudoscalar bosons in our 
model: 

A. M ~ vt- The mass of is expected to be greater than 200 GeV if m < 1 TeV. 
But if m is much larger than that, the mass of -Pj^^ will saturate to a constant 
value which is at most 0(600) GeV; the P^^ states are expected to be very heavy 
in the large m limit. The singly charged Higgs has a mass at the weak scale that is 
m-independent. The masses of neutral Higgs-like bosons are also of the weak scale. 
Being the would-be Majoron, provides a bound on uj: \uj\ > 0.016. 

B. M ~ t> ^ Vt- Here, only the mass of P^^ is expected to be at the weak scale. All 
the other scalars with the exception of (which is mostly a SM Higgs boson) will 
be too heavy to be of interest at the LHC, since their masses are controlled by uj. 



3. Neutrino phenomenology and constraints 

3.1 Two-loop neutrino masses and neutrino oscillations 

A feature of our model is that neutrino masses are generated at the two-loop level, and the 
crucial couplings are the Yukawa terms. It is well known that the Yukawa couplings of (f) to 
fermions are diagonalized by a biunitary transformation effected via Ul and C//j such that 
the charged leptons are mass eigenstates. Clearly, applying this transformation does not 
in general diagonalize 1^^. Hence we expect flavor violating couplings between families 
of RH leptons and the physical P~^^ states. ^ Thus in general, the decay modes such as 



'e'^ must occur. The coupling of P^ to fermions, on the other hand, is similar 



to SM but scaled by a factor vt/v. 

(a) 



(b) 



VaL 










1 *++ \ 








' X— 


— ' ¥i 





La la 



Figure 3: The two-loop diagram for the neutrino mass in: (a) the mass eigenbasis and (b) the 
weak eigenbasis. 

The active neutrino mass matrix can now be calculated. The leading contribution is 
given by the two-loop Feynman diagram depicted in Fig. ^. After a standard but lengthy 



^In the following we assume that the charged leptons are in the mass basis, and so Y' = UrYUr. For 
notational simplicity we will drop the prime henceforth. 
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calculation we find 



{mu)ab = -^g^mambVTYabsm{26) [l{M^r,Mp^,ma,mb) 

where a,b = e, /i, r. The integral / is given by 
I{M^,Mj,^,ml,ml) = 



,ma,mb) , (3.1) 



(2vr) 



(27r)4 P - m2 fc2 - - q^ - {k - qf - M|.. 



(3.2) 



The integral can be evaluated analytically as in [^. Note that there is a generalized GIM 
mechanism at work here. This can be seen clearly in the limit Mp^ ^ > Mw llCj: 

1 1 



I {Mi,, Mi, 0,0) 



(47r)4 M|,^ 



log' 



Ml 



(3.3) 



We see that not only is the neutrino mass two-loop suppressed, there is also a helicity 
suppression from the charged lepton masses, whose origin can be clearly seen from Fig. ^(b). 
It is clear that the internal lepton lines must have mass insertions since ^'++ only couples to 
RH leptons. As a result, {my)ee will be vanishingly small. This has important consequences 
in 0v(3(3 decays as well as the choice of signatures for the detection of these scalars at the 
LHC. 

Explicitly, the neutrino mass matrix is given by 



= f{Mp„Mp^) X 



f{Mp„Mp,) X 



/ mlYee niem^Ye^ memrYer\ 

meVn^Ye^ ml Y^^ nirm^Y^r 

ymeUlrYer nirm^Y^r mlYrr J 
/2.6 X lO-'^Yee 5.4 X 10-^ Yefj, 9.1 X 10~^Ye, 



5.4 X 10-5 Fe^. 1.1 X 10-2 

\9.lX 10-^ Yer 0.19 y^^ 



0.19 y^, 
3.i7y^^ 



(3.4) 



where 



f{Mp„Mp, 



V2g%T sin(2(5) 
128^ 



log' 



\MpJ 



log- 



■(Mw\ 
\MpJ 



(3.5) 



and f = f X (IGeV^) gives a qualitatively estimate of the overall scale of active neutrino 
masses. Now for normal hierarchy, the neutrino mass matrix has the following structure 



(3.6) 

6[), we see there is a qualitatively 




where e, e' and r] <^ 1. Comparing Eq. 
agreement. 



to Eq. 
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/(eV) 




Figure 4: The scale f{Mp-^, Mp^) as a function of Mp^ and the mass difference AAf ~ Mp^ — Mp^ 
for sin 2(5 = 0.5 and wj- = 4 GeV. 



/(eV) 
. 




Figure 5: The scale f{Mp^ , Mp^ ) as a function of Mp-^ and sin 2(5 for Mp^ — 1 TeV and vt — ^ GeV. 



We plot the behavior of f{Mp^,Mp^) as a function of Mp^ and AM = Mp^ - Mp^ 
in Fig. The range of parameters used is apphcable to Case A. We see that overall, the 
neutrino mass increases as the mass difference of the two doubly charged scalar increases. 
In Fig. ^ we plot the behavior of f{Mp^,Mp^) for parameter range applicable to Case B. 
In both cases we expect neutrino masses to be in the sub-eV range. 

We proceed next to examine how the neutrino oscillation data can constrains our 
model. Neutrino oscillations depend on the difference of mass-squared, hence we will focus 
on accordingly. Since the eigenvalues of are in general complex, and it is customary 
to separate out a phase matrix, we write 



m,. 



PMNS 



u 







1712 



U VpMNS ■ 



u 



rrig 




(3.7) 



where Vpmns is the neutrino mixing matrix |11] and in standard notation is the same as 
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the quark mixing matrix Q, whereas a is the Dirac phase. For normal hierarchy, we have 
m\ ~ 0, m| ~ Ami,, and ~ Am^^^. Oscillation experiments currently place limits on 
the following relevant parameters 



7.1 X 10"^ < Am^ < 8.9 x 10~^ (eV^) , 



0.164 < sin2 6'i2 < 0.494, 



1.4 X 10"^ < 



atm 



< 3.3 X IQ-^ (eV^) , 0.22 < sin^ 6^23 < 0.85 , 
sin^ 26*13 = ± 0.04 . 



(3.8) 



Using Eqs. (3.4) and ( |3.7| ) and the oscillation data we can get six constraints on the elements 
of m?,. Prom the first row of m?, we obtain 



fY^^ < 1.32 X lO^eV^ 



f YerYrr < 9.0 X 10"^ eV^ . (3. 



Since the phases involved are unknown, we do not get lower bounds for these quantities. 
As will be seen below, the bounds given in Eq. ( |3.9| ) are very loose compared to that 
found from rare muon and r decays. The remaining three constraints are more stringent 
relatively, and they came by demanding a good fit to normal hierarchy: 



MM 



300y^2^) < 2.7eV2 



f\Y^,, + 285y,,)y^, < 2.4 X 10^1 eV^ , 
f{Yl^^ + 278^,2^) < 2.9 X 10-2 eV^ . 



(3.10) 
(3.11) 
(3.12) 



Using Eqs. ( ^.11 ) and ( ^.12 ), we plot in Fig. ^ the allowed parameter space for 1^,- and 
Yrr, with / set to 0.5 eV. ^ 




0.1 0.2 0.3 0.4 0.5 0.5 0.7 



Figure 6: Parameter space allowed for Ytt and Y^t for / = 0.5 eV. The hyperbolic curve is the 



upper limit set by Eq. ( 3.11 ). The shaded region is that allowed by both Eq. ( 3.11 ) and Eq. ( 3.12 ) 



*In plotting Fig. ^ we have used the fact that Y^^Y^^ <C 0.24, which we show below. 
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3.2 Rare muon and r decays 

The doubly charged Higgs bosons lead to many lepton number violating processes. Since 
no such signals were found in current experiments, they lead to very strong constraints on 
the Yukawa couplings. In the following, we work out these constraints. 

A. Muonium anti-muonium conversion 

The effective Hamiltonian is given by the exchange at tree level: 

Y Y 

Hmm = 2m2^ P-l^f^R P-li^eR + h.c. , (3.13) 

where M is the reduced mass of the pair of doubly charged Higgs given by 

1 sin^ 6 cos^ 5 ,„ ^ 

+ -Trr- ■ (3-14) 



The current experimental limit B gives 



YeeY^^ < 2.0 X 10"^ (M__/100GeV)2 . (3.15) 



B. EfTective e+e l^l , I = e,fi,T, contact interactions 

The effective Hamiltonian for Bhabha scattering is 



Y'^ _ 

eRj^en enjf.eR . (3.16) 



The bounds are 



Yl < 1.8 X 10^^ (M__/100GeV)2 , 
Y^^ < 2.4 X 10^3 (A/__/100GeV)2 , 

Y^^ < 2.4 X 10^3 (M__/100GeV)2 . (3.17) 

C. Rare ^ ^ 3e decays and its r counterparts 

These decays can all be induced at the tree level and thus provide the most stringent 
limits on the Yukawa couplings. For fj, — > 3e, the branching ratio is given by 

«.„.3.).(1^)^(^)\ (3.18, 

with similar equations for r decays. The constraints impose by the data is given by 

Yef^Yee < 6.6 X 10"^ (M__ /lOO GeV)^ 
YerYee < 3.0 X 10"^ (M__ /lOO GeV)^ 
YerYf^f, < 3.0 X IQ-^ (M__/100GeV)2 
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Yf,rYf,f, < 2.9 X 10"'* (M__/100GeV)^ 
Yf^rYee < 2.9 X 10""^ (M__ /lOO GcV)^ 



(3.19) 
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D. Radiative flavor violating charged leptonic decays 

We consider here rare radiative decays of — > 67 and r /u(e)7- The fact that 
they are not seen to a very high precision make them of paramount importance for 
probing the physics of lepton flavor violation. The branching for ^ ^ 67 is calculated 
in ||T^ using an effective theory approach: 

M^^^eJ) = ^ (3.20) 



37rG2 

^ l=e,fj.,T 



with obvious substitutions for r decays. The limits are given by 

Y^Yii^Yie < 1-5 X 10-5(M__/100GeV)2, 
I 

Y^^irYie < 1.4 X 10-3(M__/100GeV)2, 



< 1.1 X 10"^ (M__/100GeV)' 



(3.21) 



Comparing the sets of constraints we find that Eq. ( |3.19| ) gives the strongest limits. 
Although the limits from contact interactions, Eq. ( 3.17| ), are less stringent, they are useful 
nonetheless as they constrain individual couplings. We illustrate in Fig. |^ how Y^r and Y^e 
are restricted by both contact and rare decay experiments for a chosen value of the reduced 
mass M__ = 400 GeV. 




0.05 



. 1 



0.2 



0.3 



0.4 



Figure 7: Parameter space allowed for Yer and Yee (shaded). The straight lines are limits from 
the contact interactions, the hyperbolic curve the upper limit from t 3e decays. 

It is also interesting to compare the limits from the rare decays with that from neutrino 
oscillations. We plot in Fig. |8| the parameter space allowed in this case for Y^r and Y^^, 
with / = 0.5 eV and M = 400 GeV. Other comparisons are less instructive. 

In summary, we have from contact interactions the upper limits 



<0.17, 



< 0.2 , 



Y^r < 0.2 . 



(3.22) 
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10 K 



1 . 5 



. 5 




3.5 



Figure 8: Parameter space allowed for for Y^^ and Y^^ (shaded). The hyperbolic curve is the 
upper limit from t — s- 3/^ decays, the ellipse the limit from fitting neutrino data. 



and from neutrino data (see Fig. ^ and |^ 



Yrr < 0.02 . 



(3.23) 



The values M = 400 GeV and / = 0.5 eV are used throughout in obtaining these limits. 

They are consistent with and typical of what can be expected for our model. 

3.3 Ovf3P decays of nuclei 

In our model, Oi^/3/5 decays of nuclei are induced by the exchanges of virtual Pi 2 bosons 
and Majorana neutrino as depicted in Fig. ^. The quark level amplitude due to neutrino 

(a) (b) 



w- 



J,2 




en 



Figure 9: Oi^/3/? decays via exchange of: (a) doubly charged Higgs and (b) light Majorana neutrinos. 



exchanges is given by 



(3.24) 



where < p > is the average momentum of the light neutrino exchanged. For notational 
simplicity we will drop the subscript v for the mass matrix. Typically, < p >~ 0.1 GeV 
which reflects the long range nature of light particle exchanges. Note that there is a 



^In our order of magnitude estimation, we have ignored spinor and kinematic factors, as well as factors 
from nuclear physics. 
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cancellation between the contributions from Pi and P2, which is characteristic of a two 
level system. 

The doubly charged Higgs exchange amplitude is given by 



A 



(7^ YeeVT sin 2(5 / 1 



1 



w 



(3.25) 



where niee is given by Eq. (^.4|). We estimate that A,y/Ap-- < 10 ^. The smallness of 
this ratio is due to the fact that in our model, rriee is suppressed not only by a two-loop 
factor, it is also suppressed by the electron mass factor (mg/Miy)^ coming from the doubly 
charged scalar coupling. We conclude that if seen, OvPf] decays of nuclei will be due to the 
existence of doubly charged Higgs at the weak scale. This can be tested at the LHC, and 
is the subject of next section. Since there is no conclusive evidence for these decays, we 
use it to set a limit of Y^e- The result is displayed in Fig. It can be seen that the limit 
here is comparable to that from the contact interactions. 

Y 

p.p. 



0.4 



0.3 
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400 
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Figure 10: Upper limit on as a function of A/p^ for sin 25 = 0.5 and Mp^ = 1 TeV. 



4. Doubly Charged Higgs at the LHC 
4.1 Production of the doubly charged Higgs 

A central ingredient in our neutrino mass generation is the two doubly charged Higgs, -Pf^2^. 
We have argued that if not both, at least one of the doubly charged Higgs is well within 
reach of the LHC. Without loss of generality, we will take Pi to be this (lighter) state, and 
focus on its production and decay below. 

Now the doubly charged Higgs, P^^, have no direct couplings to the quarks, which is 
characteristic of models with doubly charged scalars. However, they do couple to the SM 
gauge bosons. Thus at the LHC, Pj^^ will be produced predominantly via the WW fusion 
processes, as illustrated in Fig. ^(a) (the d-quark is replaced by the u-quark for Pi"*"*"), and 
the Drell-Yan (DY) annihilation processes, 

qq^^^Z* ^P++P{- {q = u,d). (4.1) 
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The relevant gauge-scalar couplings are given by 



P^+P^- ■.i2eA^ d^P++P{- + h.c. 



- ■ [{l-2sl,)cj-2sl,sj]Z,d,P++P,- 



h.c. 



cw 



(4.2) 



where g = e/ sinOyy, sw = sinOw, cs = cos 5 and ss = sin 5. 

We calculate the production cross sections numerically using the CalcHEP package |13|, 
which allows an easy implementation of our model. The calculations are leading order, and 
are done in unitary gauge. The cross sections are calculated for pp collisions in the center- 
of-mass (CM) frame with energy ^/s = 14 TeV, and CTEQ6M |14] parton distributions 
functions are used to fold in the cross section for the hard partonic production processes. 

From Eq. (O) we see that the only model parameters the production cross section 
explicitly depend on are vt and the mixing angle 6. The dependence on all other model pa- 
rameters are implicit through the dependence on the -Pf*^^ mass, Mp^, as given in Eq. (p.^). 
The choice we made in our calculations is thus to set 



M = VT, Kq, = -K2, m = 2u = 492.442 GeV, A = -A^^ = p/2 = 1 



(4.3) 



and vary Mpj by varying K2. Note that with this choice, ss is independent of K2, and thus 
Mpj. For the SM parameters, we take v = 246.221 GeV, e = 0.3133 and sw = 0.4723. 

a(fb) 




200 



300 



400 



Figure 11: Production cross sections for with vt ~ 4, GeV and ss = 0.12. The solid line 

traces the results from the W fusion processes, the dashed line the Droll- Yan annihilation processes. 
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We plot in Fig. 11 the production cross section for P^^ from WW fusion and Drell-Yan 
pair annihilation as a function of Mp^ . Note that QCD corrections are expected to increase 
the DY production cross section by a factor of about 1.25 at next-to- leading order |15|. 
The CalcHEP results are checked with that calculated from Pythia [16|, and are found to 
be consistent. 



^For couplings with P2, 



cs 



cs- 
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Except for the mixing angle factors cs and s^, the coupling of to the weak gauge 
bosons are the same as that of to the Wl, Z\ bosons in the left-right symmetric 



model 1 17], while the couplings to photon is the same in both models. Thus, the results for 
the cross section should agree once the scaling factors are taken into account. Comparing 
the results of our model to that of Refs. we find good agreements. 

We see from Fig. |l^ that DY pair production is dominant compared to the single 
production channel via WW fusion for the whole range of Pf'^ masses in our model. This 
is a robust feature in the production of the doubly charged Higgs at hadron colliders in 
many models. This is because the virtual photon exchange is the same for all models 
and only the Z exchange term is model dependent, which is subdominant. The relative 
magnitude of production via WW fusion compared to the DY pair production can thus be 
used to distinguish between the different models containing the doubly charged Higgs. 

4.2 The decay of P^^ 

At the leading order, there are six decay channels possible for P^^: 

(1) P±± ^ (a,6 = e,/i,r), 

(2) P^^ ^W^W^ , 

(3) Pf^ ^P^W^ , 

(4) P^^ ^ P^P^ , 

(5) P^^ ^ W^W^X^ , xo = ro, h^, po 

(6) ± ^ P^P^X^ . (4.4) 

Kinematically, mode (4) and (6) are not allowed in our model, while the availability of the 
rest depends on the value of the scalar boson masses. The coupling for mode (2) has been 
given in Eq. ([4.2|), for mode (1), (3) and (5) the couplings are given by 

Pl^l^R^^R Yab SsPr'^hR + h.C. 

Pf^W;fP^ : ig cs W~ [d^P^+P" - P^+d^P^] + h.c. 

Pt^W;fW;^X^ : ^ cs cxPt^W;fW;fx'^ + h.c. (4.5) 
v2 

where cx = i, cosi?, sini? for X^ = T^,h^,P^ respectively. 

The leptonic decays are kinematically the most favorable mode of decay for the doubly 
charge Higgs, -Pf*^^, and this is a universal feature in any model that contains them. The 
leptonic decay width has a very simple form given by 

r(44) = (l + '^a6)^s|Mp, (no sum). (4.6) 

Note that in our model, the final state charged leptons are right-handed. Hence in principle, 
helicity measurements can be used to distinguish between our model and those whose 
doubly charged Higgs coupling only to left-handed leptons (see e.g. |2l|, |2^, p4|). 

From the discussion above, it is not unreasonable to take Yab ~ 0.1 (except for Yrr, 
which is constrained to be ten times smaller). Thus, provided the mixing angle between 
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the doubly charged Higgs bosons, 6, is not too smaU, one can expect spectacular signals 
from like-sign dileptons such as e^, er and fj,T, that are directly produced from i^j^^- 
The l^^H^^ channel opens up once Mp^ > 2M\y. Its decay width is 



levrMp, 



4M^ 



(4.7) 



which is proportional to v^. Given that vt is small, and Mp^ is of order the electroweak 
scale, whether the leptonic or W^W^ mode dominates depends on the value of the mixing 
angle 5. 

5M (GeV) 
150 




600 



Mp, (GeV) 



Figure 12: The mass difference 6M = Mp^ — Mp as a function of Mp^ for = 4 GeV. The short 
dashed line marks where 5M = Mw, which opens up the decay channel. 



The W^P^ mode opens up once the mass difference, 5M = Mp^ — Mp, where Mp 



is the mass of the singly charged Higgs boson, is greater than 2M\y. From Eq. (2.7), we 
see that if Mp is to be of interest at the LHC, w needs to be order unity (the validity of 
perturbation theory constrains \k2\ to be less than 47r). This implies that M ^ vt, and in 
Fig. |l^ we plot 5M as a function of Mp^ for M = vp = ^ GeV. We vary Mp^ the same 



way as detailed in Sec. 4.1, with the same choice of parameters, Eq. ([4.3|). For this choice, 
the decay channel opens up once Mp^ > 407 GeV. 
The decay width for the W^P^ mode is given by 



T{W^P^) 



g'cjMl 



WttM^ 



I Mk 

'Ml' 



Ml 



(4.8) 



where A(x, y, z) = + + — Ixy — 2xz — 2yz. With a dependence of and no 
suppression coming from factors of vp, the W^P^ mode is expected to dominate over 
the W^W^ mode once it becomes available. Note that the singly charged Higgs decays 
primarily into rather than fermion pairs, which is the dominant decay mode for the 
usual charged Higgs originating from SM doublets. 

In general, the three-body decay mode W^W^X^ , = T^,h^,P^, is expected to 
be relatively suppressed by the three-body phase space when compared to the two-body 
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modes. However, since the couplings do not depend on vt, if ss <C 1, it may become 
important compared to the lepton and W^^M^^ modes. 

The masses of the neutral pseudoscalar and scalar bosons, nixo , are given by Eqs. ( 2.1^ ) 
and ( 2.15 ). Here too, uj ~ 0(1) is required to keep the neutral bosons at electroweak scale. 
Note that once to is fixed, Mxo and Mp^ are governed by independent sets of parameters. 
Below, we will choose = as the representative three-body mode, which has the least 
and simplest model dependence compared to the neutral scalar cases. The difference in the 
decay widths would come in only as ratios of masses, and from the neutral mixing angle 
factors, and s^. 
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Figure 13: Decay widths of -Pj^^ as a function of Mp^^ with = 4 GeV and the mixing angle 
fixed at sin (5 = 0.12. The leptonic mode is represented by the solid lines at different values of 
l^^abP, the W^W^ mode the long dashed line, and the W^P^ mode the short dashed line. The 
three-body W^W^T^ mode is labelled on the plot. 



For the choice M = = 4 GeV, M^o = 146.5 GeV, and the W^W^T^^ mode opens 
up once Afp^ > 307.31 GeV. Since there is no known closed form for the phase space of 
general three-body decays, we calculate the three-body decay widths numerically using 
CalcHEP as a function of Mp^. We display the result and compared with all the allowed 



two-body modes listed in Eq. (4^) in Fig. 13. We work in the small mixing scenario where 
Ss = 0.12, since for ~ 1, the results are qualitatively that of the left-right symmetric 
models 0. 

There are several noteworthy features. First, despite the phase space advantage, the 
suppression of the lepton mode by the s| factor is such that the VF^T^^ mode is never 
negligible for \Yab\'^ ~ 1, and dominating for <C 1 even for small values of Mpj. Next 

as discussed above, the W^P^ mode dominates over both the lepton and T^^Vl^^ modes, 
and does so almost at once after it becomes kinematically available. Moreover, even if one 
were to push jYabP to the perturbative limit of Air, the W^P^ mode would still become 
dominant soon above its threshold. Lastly, and perhaps most interestingly, the three-body 
W^W^T^ mode dominates over the lepton mode when l^afeP ^ 1 (suppression by ss), 
and is significant relative to the mode and even overtakes it at large values of Mp^ 

(suppression by vt)- 
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These are features that are not shared by other doubly charged Higgs models, and will 
serve to distinguish our model from them at the LHC. 

5. Conclusion 

We have given a detail study of the phenomenology of a model where Major ana neutrino 
mass is induced radiatively at two-loop. This is made possible by extending the Higgs 

sector of the SM, and the new structure involves a triplet and a doubly charged singlet 
field. Lepton number is violated only in the scalar sector by a hard and a soft term. When 
these two terms are set to zero, the model conserves lepton number, and is technically 
natural. 

Since the active neutrinos masses are two-loop effects due to new physics at the TeV 
scale, the model is testable at the LHC. In our study, we have focused on the production 
and decays of the doubly charged Higgs bosons. We find it is possible in the enlarged 
Higgs sector parameter space to have always at least one doubly charged Higgs with mass 
< 600 GeV. Furthermore, in much of the parameter space, two doubly charge Higgs can 
be expected to have masses within reach of the LHC. We find that Drell-Yan mechanism 
is the dominant production mechanism, with single production from WW fusion coming 
into play only for masses greater than 600 GeV. This is common to most doubly charged 
Higgs models. Wc find cross sections ranging from 0(1 — 70) fb without including QCD 
effects, which typically cause an enhancement by a factor of 1.25. 

The decays of the doubly charged Higgs are potentially the best way of differentiating 
between the different doubly charged Higgs models. We find that with reasonable val- 
ues of the Yukawa couplings, the leptonic decays are not always dominant despite being 
kinematically favorable. Nevertheless, decay modes producing directly same sign dilepton 
pairs such as er, /it, and e/i, are distinctive and unmistakable. Moreover, since the doubly 
charged Higgs only couples to RH Icptons in our model, their decays into a pair of tauons 
offer opportunities for helicity measurements. Another interesting feature of the doubly 
charged Higgs decays is that three-body modes W^W^X^ involving the neutral scalar and 
pseudoscalar bosons can become as important as the two-body modes, particularly when 
the mixing between the doubly charged Higgs is small. 

The model we studied here can only accommodate the normal hierarchy for the active 
neutrino mass matrix. It predicts that the first element to be very small. Moreover, Oz^/3/5 
decays of nuclei can still occur via the exchange of doubly charged Higgs. This give rise to 
the interesting situation whereby if such decays are seen, the doubly charged Higgs will also 
be seen at the LHC, and vice versa. If the neutrino mass hierarchy turns out to be other 
than the normal hierarchy, neutrino masses are not generated by the model we studied. 
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